

clear all
set more off
capture log close
set mat 5000

cd "C:\Users\rujia\Dropbox\coalmine\restat"

use  CoalMine.dta, clear
drop if native==.
drop if leader_id==.


gen ln_output=log(prod0)

gen prod0_normalized=prod0/10000
label var prod0_normalized "Output"


label var ln_output "ln Output"


bysort prov_id: egen pre_deathrate=mean(deathrate) if year<1998
bysort prov_id: egen m_pre_deathrate=mean(pre_deathrate) 
gen pre_year=m_pre_deathrate*year
gen decentralized=(year>=1998 & year<=2000)

tabu prov_id, gen (pdummy)
forvalues i=1(1)22{
gen prov_year`i'=pdummy`i'*year
}

gen period=1 if year<1998
replace period=2 if year>=1998 & year<=2000
replace period=3 if year>2000 & year<=2005

egen clusterg=group(prov_id period)

gen decentralized_native=decentralized*native

gen ln_wage=log(indwage)
gen tenure=year-startyear+1
gen ln_gdpper=log(gdpper)
gen ln_electricity=log(eleconsum)
gen ln_disbeijing=log(1+dis_beijing)


local controls "ln_output ln_wage tenure age ln_gdpper ln_electricity gov_native s_native ln_disbeijing"


foreach var in ln_output ln_wage tenure age ln_gdpper ln_electricity gov_native s_native ln_disbeijing{
gen decentralized_`var'=decentralized*`var'
}






gen ln_death=log(death)

drop if year>2005
tabu year, gen(ydummy)


forvalues i= 1(1)11{
gen native_year`i'=ydummy`i'*native
}
label var native_year1 "Native * Year 1995"
label var native_year2 "Native * Year 1996"
label var native_year3 "Native * Year 1997"
label var native_year4 "Native * Year 1998"
label var native_year5 "Native * Year 1999"
label var native_year6 "Native * Year 2000"
label var native_year7 "Native * Year 2001"
label var native_year8 "Native * Year 2002"
label var native_year9 "Native * Year 2003"
label var native_year10 "Native * Year 2004"
label var native_year11 "Native * Year 2005"


drop decentralized_native



#delimit ;

eststo clear;
eststo:  reg deathrate native_year1 native_year2 native_year4-native_year11, cluster(clusterg);

estadd scalar nobs=e(N);
estadd scalar rsq=e(r2);




eststo:  xi:reg deathrate native_year1 native_year2 native_year4-native_year11 i.prov_id i.year, cluster(clusterg);
estadd local pfe "Y";
estadd local yfe "Y";

estadd scalar nobs=e(N);
estadd scalar rsq=e(r2);


eststo:  xi:reg deathrate native_year1 native_year2 native_year4-native_year11 `controls' i.prov_id i.year , cluster(clusterg);
estadd local pfe "Y";
estadd local yfe "Y";
estadd local controls "Y";

estadd scalar nobs=e(N);
estadd scalar rsq=e(r2);

eststo:  xi: reg deathrate native_year1 native_year2 native_year4-native_year11 `controls' decentralized_* pdummy* ydummy*, cluster(clusterg);

estadd local pfe "Y";
estadd local yfe "Y";
estadd local controls "Y";
estadd local controls_d "Y";

estadd scalar nobs=e(N);
estadd scalar rsq=e(r2);





esttab using draft/Dynamic.tex, replace
	width(\hsize) 
	nomtitles
	b(3)
	se(3)
	noconstant
	star(* 0.10 ** 0.05 *** 0.01)
	keep(
native_year1 
native_year2 
native_year4
native_year5
native_year6
native_year7
native_year8
native_year9
native_year10
native_year11
		)
	label
	stats(
		pfe 
		yfe
		
		controls
		controls_d

		nobs
		rsq
		, 
		labels(
			"Province FE" 
			"Year FE"
	
			"Controls"
			"Decentralization * Controls"

			"\# observations"
			"R-Squared"
			) 
		fmt(0 0 0 0 0   2)
		layout(@)
		)
	nonotes
	addnotes(
	)
	;

